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We measure branching fractions and integrated rate asymmetries for the rare decays B — » 
K^*'t^(.~, where l'^l~ is either e^e^ or fi'^ /i^ , using a sample of 384 million BB events col- 
lected with the BABAR detector at the PEP-II e+e^ collider. We find no evidence for direct CP 
or lepton-flavor asymmetries. However, for dilepton masses below the J/^ resonance, we find ev- 
idence for unexpectedly large isospin asymmetries in both B — > Kt^ i~ and B — > K*£'^£~ which 
differ respectively by 3.2a and 2.7a, including systematic uncertainties, from the Standard Model 
expectations. 
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The decays B ^ K'^*H+(.-, where £+£" is either e+e" 
or fi^ fi^ , arise from flavor-changing neutral current pro- 
cesses that are forbidden at tree level in the Standard 
Model (SM). The lowest-order SM processes contributing 
to these decays are a W~^W~ box diagram, and the ra- 
diative photon and electroweak Z penguin diagrams [l|. 
Their amplitudes are expressed in terms of hadronic form 
factors and effective Wilson coefficients C^ , C^ and 
C^Q , representing the electromagnetic penguin diagram, 
and the vector part and the axial-vector part of the Z 
penguin and W'^W~ box diagrams, respectively [3|- New 
physics contributions may enter the penguin and box di- 
agrams at the same order as the SM diagrams, modifying 
the Wilson coefficients from their SM expectations [3|. 

We report results herein on exclusive branching frac- 
tions, direct CP asymmetries, the ratio of rates to di- 
muon and di-electron final states, and isospin asymme- 
tries, measured in two regions of dilepton mass squared 
chosen to exclude the region of the J/ij) resonance: a low 
q^ region 0.1 < q^ = mj^ < 7.02 GeV^/c'' and a high q'^ 
region q^ > 10.24 GeV^c^. We also present results for 
the two regions combined. The ipi'^S) resonance is re- 
moved from the high q^ region by vetoing events with 
12.96 < q^ < 14.06 GeV^/c"*. For K*e+e- final states, 
we also report results in extended low and extended com- 
bined q^ regions including events g^ < 0.1 GeV^/c*, where 
there is an enhanced coupling to the photonic penguin 
amplitude unique to this mode. Recent BaBAR results 
on angular observables using the same dataset and simi- 
lar event selection as is used here are reported in [l| . 

The B — > K£~^i~ branching fraction is predicted to 
be (0.35 ± 0.12) x 10-^ while B ^ K*i+e- for q"^ > 
0.1 GeV^c^ is expected to be roughly three times larger 
at (1.19 ± 0.39) X 10"*^ 0. The ~ 30% uncertainties 
are due to lack of knowledge about the form factors that 
model the hadronic effects in the B ^ K and B -^ K* 
transitions. Thus, measurements of decay rates to exclu- 
sive final states are less suited to searches for new physics 
than rate asymmetries, where many theory uncertainties 
cancel [J]. 

The direct CP asymmetry 
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is expected to be 0{lO^^) in the SM, but new physics 
at the electroweak scale could produce a significant en- 
hancement [5|. 

The ratio of rates to di-muon and di-electron final 
states 
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is unity in the SM to within a few percent [6|. In two- 
Higgs-doublet models, including supersymmetry, these 
ratios are sensitive to the presence of a neutral Higgs 
boson, which might, at large tan/3, increase Rj^m by 
~ 10% |7[. In the region q"^ < (2m^)^, where only the 
e+e~ modes are allowed, there is a large enhancement 
of i? — > K*e^e~ due to a l/q"^ scaling of the photon 
penguin. The expected SM value of Rk* including this 
region is 0.75 [6|, and we fit the K* dataset over the ex- 
tended combined and extended low q^ regions in order to 
test this prediction. 

The CP-averaged isospin asymmetry 
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where r = rn/r+ = 1/(1.07 ± 0.01) is the ratio of the B° 
and i?+ lifetimes l8[ , has a SM expectation of -1-6— 13% as 
q^ -^ GcY^/c^ |9[. This is consistent with the measured 
asymmetry of 3±3% in B -^ K*^ [8|. A calculation 
of the predicted K*'^ and K*'^ rates integrated over the 
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low q^ region gives Af 
the high q^ region, contributions from charmonium states 
may provide an additional source of isospin asymmetry, 
although the measured asymmetry in J/ipK^*' is at most 
a few percent (8( . 

We use a data sample of 384 million BB pairs collected 
at the T(4iS') resonance with the BABAR detector [12| at 
the PEP-II asymmetric-energy e^e~ collider at SLAG. 
Our selection of charged and neutral particle candidates, 
as well as reconstruction of 7r°, K^ and K* candidates, 
are described at [l| . We reconstruct signal events in ten 
separate final states containing an e+e~ or /i^/i~ pair, 
and a K^[-^ tt+tt"), K^ , or K*{%Q2) candidate with an 
invariant mass 0.82 < M{K'n:) < 0.97GeV/c2. We recon- 
struct K*^ candidates in the final state K~^Tr~ , and K*~^ 
candidates in the final states /\+7r° and K'^tt'^ (charge 
conjugation is implied throughout except as explicitly 
noted). We also study final states K^*^h^fj.^, where ft, is a 
track with no particle identification requirement applied, 
to characterize backgrounds from hadrons misidcntificd 
as muons. 

B -^ K^*^£'^£~ decays are reconstructed using the 
kinematic variables 7ties ~ y^s/4 — p*^ and AE = 
Eg — \/s/2, where pg and Eg are the B momentum 
and energy in the T{AS) ccnter-of-mass (CM) frame, 
and Y^ is the total GM energy. We define a fit region 
mss > 5.2GeV/c2, with -0.07 < AE < 0.04 (-0.04 < 
AE < 0.04) GeV for e+e" (a'^M^) final states in the low 
and extended low g^ region, and —0.08 < AE < 0.05 
(-0.05 <AE< 0.05) GeV for high q^. 

The main backgrounds arise from random combina- 
tions of leptons from semileptonic B and D decays, which 
are suppressed through the use of neural networks (NN) 



whose construction is described in detail in [l| . For each 
of the ten final states we use separate NN optimized to 
suppress either continuum or BB backgrounds in the low, 
extended low or high q^ regions. We use simulated sam- 
ples of signal and background events in the construction 
of the NN, and assume rates consistent with accepted 
values [8|]. 

There is a further background contribution from B — > 
D{-^ K'-*'>'!r)TT decays, where both pions arc misidcnti- 
fied as leptons. The pion misidentification rates are 2-3% 
for muons and <0.1% for electrons, so this background 
is only significant in the /x"'"/i~ final states. We veto 
these events by assigning the pion mass to a muon candi- 
date, and requiring the invariant mass of the hypothetical 
K^*^TT system to be outside the range 1.84-1.90 GeV/c^. 
After all the above selections have been applied, the fi- 
nal reconstruction efficiency for signal events varies from 



3.5% for K+t:^ 



for K 



-+-, 



in the high q^ region. 



for the combined q region, to 22% 



We perform unbinncd maximum likelihood fits to toes 
distributions to obtain signal and background yields. We 
use an ARGUS shape [l3| to describe the combinatorial 
background, allowing the shape parameter to fioat in the 
fits. For the signal, wc use a fixed Gaussian shape unique 
to each final state, with mean and width determined from 
fits to the analogous final states in the vetoed J/V-'A'*-*' 
events. We account for a small residual contribution 
from misidentified hadrons by constructing a probability 
density function (pdf) using K^*^h^iJ,^ events weighted 
by the probability for the h^ to be misidentified as a 
muon. We also account for background events that peak 
in the niES signal region, arising from charmonium events 
that escape the veto, and for contributions from misre- 
constructed signal events. We test our fits in each fi- 
nal state using the large samples of vetoed J/t/iAT^*' and 
^{2S)K^*'^ events, and find that all the branching frac- 
tions are in good agreement with accepted values [14| . 

Wc perform simultaneous fits for A^p , Rxi*) and Af 
summed over all the signal modes that contribute to the 
particular measurement. 

We estimate the statistical significance of our fits by 
generating ensembles of 1000 datasets for each of the ten 
final states in each q^ region of interest, and fitting each 
dataset with the full fit model described above. These 
tests also confirm the unbiased nature and proper error 
scaling of our fit methodology. 

For the total B -^ Kl+i' and B -^ K*e+£- branch- 
ing fractions averaged assuming isospin and lepton-flavor 
symmetry, we measure (0.394to ^gg ±0.020) x 10^^ and 
(l.lltai^ ±0.07) X 10"^ respectively, where the first 
uncertainty is statistical and the second is systematic. 
Complete branching fraction results in all final states and 
q^ regions, along with the statistical significance of each 
measurement and frequcntist upper limits for measure- 
ments with < 4(T statistical significance, are available on- 



line [15[ . All results are in good agreement with previous 
measurements [8|. 

Table |T] summarizes the results for A^p . In the fits 
to the separate B and B datasets in charge-conjugate 
final states, we assume a common background ARGUS 
shape parameter. Our final results are consistent with 
the SM expectation of negligible direct CP asymmetry. 
Table |lT] shows the results for Rk and Rk* , which are 
also consistent with the SM expectations. 

Table IIIII shows the results for the isospin asymme- 
try Af . We directly fit the data for Af taking 
into account the differing lifetimes of S" and B'^ . Fig- 
ure [T] shows the charged and neutral low q^ datasets with 
overlaid fit projections. We find no significant isospin 
asymmetries in the high and combined q^ regions, or for 
K*e^e^ fits in the extended regions. However, we find 
evidence for large negative asymmetries in the low q^ re- 
gion. 

We calculate the statistical significance with which a 
null isospin asymmetry hypothesis is rejected using the 
change in log likelihood v2Aln£ between the nominal 
fit to the data and a fit with Af — fixed. Figure [2] 
shows the likelihood curves obtained from the K£'^£~ and 
K*£'^i~ fits. The parabolic nature of the curves in the 
Af > — 1 region demonstrates the essentially Gaus- 
sian nature of our fit results in the physical region, and 
the right-side axis of Figure [2] shows purely statistical 
significances based on Gaussian coverage. Incorporating 
the relatively small systematic uncertainties as a scaling 
factor on the change in log likelihood, the significance in 
the low q^ region that Af is different from zero is 3.2a 
for Ke+£- and 2.7cr for K*l+£-. Wc have verified these 
confidence intervals by performing fits to ensembles of 



TABLE I: Aqp results in each relevant q^ region. The un- 
certainties are statistical and systematic, respectively. 



Mode 



combined q 



low q 



high q 



-0.18 ■ 

K''°i+i- 0.02+°;^° ± 0.02 



-0.19 -^ "-"^^ "•"='-0.39 

-0.231q 3g ± 0.02 0.17+°:m =•= 0-02 



K* 
K" 



O.Ollo;24 ±0.02 0.10lo:24±0-02 -0.18to55 ±0.04 



,+0.16 



,+0.21 



±0.01 0.09l°:2l ±0-02 



TABLE II: Rjfi*) results in each cf region. The extended 
("ext.") regions are relevant only for Rk* ■ The uncertainties 
are statistical and systematic, respectively. 



q^ Region 


Rk* Rk 


combined 


l-37l^jg ± 0.09 0.961^-34 ± 0.05 


ext. combined 


l.lOto.S ±0-07 — 


low 


1.0ll°j^ ± 0.08 0.40+°J° ± 0.02 


ext. low 


O-SetoJa ± 0.04 — 


high 


2.15+^-78 ±0.15 1-06+0 ji ±0.06 



TABLE III: Af results in each q^ region. The uncertainties 
are statistical and systematic, respectively. The last table row 
shows K*e'^e~ results for the extended regions. 

Mode combined q^ low q^ high q^ 

Ke+e- 

Ke+e- 

K*e+e- 

K*£+e- 



0.13+q;37±0.04 



-0.91"* 



±0.18 



-0.73 
-0.37 



+0.39 
-0.50 
+0.27 



0.39+°;4g ± 0.04 



±0.04 -1.41" 
±0.04 



0.69 



±0.04 0.21 



+0.32 



±0.03 



-O.OOtoll ± 0.05 
-0.20to:2o ± 0-03 



-0.26 
-0.66 



-0.85 

«-J«±0.05 

+0.19 



-0.12 



+0.18 



±0.04 -0.56 



-0.17 
+0.17 



±0.02 
±0.03 



0.41 

0.28to;3o ± 0.03 
-O.OStoJ? ± 0.05 
0.32+°Jg ± 0.03 
O.lSto'al ±0.04 



K*e+e- 



-0.27l°;^^±0.03 



-0.25t°o^l ± 0.03 




FIG. 2 

[long dash], K*£ 



Low q^ region Af fit likelihood curves. Kt^l 



-r [short dash], {K,K*)e.+r [solid]. 
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FIG. 1: Charged and neutral fit projections in the low q^ re- 
gion. Total fit [solid], conibinatoric background [long dash], 
signal [medium dash], hadronic background [short dash], 
peaking background [dots]. 



simulated datasets generated with Af = fixed, and 
we find frequentist coverage consistent with the A In £ 
calculations. The highly negative Af values for both 
K£^£~ and K*£~^£~ at low q^ suggest that this asymme- 
try may be insensitive to the hadronic final state, and so 
we sum the likelihood curves as shown in Figure [2] and 



^Q i[4±0.03. including systematics. 



obtain ^f **' = -0.64 

this is a 3.9(T difference from a null Af ' hypothesis. 

We consider systematic uncertainties associated with 
reconstruction efficiencies; hadronic background param- 
eterization in di-muon final states; peaking background 
contributions obtained from simulated events; and pos- 
sible CP ^ lepton flavor and isospin asymmetries in the 
background pdfs. We quantify the efhciency systemat- 
ics using the vetoed J/^'i^T^*^ samples. These include 



charged track, tt*^, and K^^ reconstruction, particle iden- 
tification, NN selection, and the AS and K* mass selec- 
tions. The largest contributions to the systematic uncer- 
tainties on the rates are particle identification, the char- 
acterization of the hadronic background and the signal 
JTiES pdf shape. All of these cancel at least partially in 
the rate asymmetries, and the final systematic uncertain- 
ties are small compared to the statistical ones. 

We perform several additional checks of effects that 
might cause a bias in our final results. We vary the pa- 
rameterization of the hadronic background pdfs, and of 
the random combinatorial background ARGUS shapes 
in the low q^ region, to test the robustness of the large 
Af asymmetries. We remove all the NN selections, 
and perform separate fits to the two K*^ final states, 
and observe no significant variation in the Af re- 
sults. To understand if an isospin asymmetry might be 
induced by the combinatorial background, we compare 
data and simulated background events within a larger 
region jASj < 0.25 GeV outside our Ai5 selection win- 
dow and in the 5.2 < toes < 5.27 GeV/c^ region. We find 
that the numbers of simulated and data events in this 
larger region agree well. No signal isospin asymmetry is 
found using simulated events within the fit region. 

In summary, we have measured branching fractions, 
and studied direct CP violation, ratios of rates to di- 
muon and di-electron final states, and isospin asymme- 
tries in the rare decays B — > K^*'>£^i^. Our branch- 
ing fraction results agree with both SM predictions and 
previous measurements. Our results for the direct CP 
asymmetries and Icpton-flavor rate ratios are in good 
agreement with their respective SM predictions of zero 
and one. The isospin asymmetries in the high and com- 
bined q- regions arc consistent with zero, but in the 
low g2 region in both B -^ K£+£- and B -^ K*£+£- 
we measure large negative asymmetries that are each 
about 3cr different from zero, including systematic un- 
certainties. Combining these results, we obtain Af — 
—0.64^0 :[4 ± 0.03, with a 3.9(t difference (including sys- 



tematics) from A^ = 0. Such large negative asymnie- article. For more information on EPAPS, 

tries are unexpected in the SM, which predicts essentially ,http://www.aip.org/pubservs/epaps.html, 

no isospin asymmetry integrated over our low q^ region 
and, as q^ — > 0, an asymmetry of ^ +10%, opposite in 
sign to our observation in the low q^ region. 
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Table [J shows the total and partial branching fractions in i? — > K^*^i'^£~ modes, along with statistical and system- 
atic errors, and the statistical significance of each measurement. To calculate significance, the data are refit with a null 
signal hypothesis and the change in ln£ with the nominal fit is used to determine the significance. For measurements 
with a statistical significance less than 4.0tT, wc fit ensembles of toy datasets derived from our observations in the data 
to compute frequentist upper limits at the 90% confidence level, in which only fit results which give a physical signal 
yield are used. The combined g^ results are scaled to account for the regions of the total dilepton mass distribution 
which are vetoed and not included in the fits. 

TABLE I: Branching fractions by mode and q^ region with statistical and systematic errors, respectively. The statistical 
significance of each measurement is given parenthetically. 90% confidence level upper limits are provided for measurements 
with < 4cr statistical significance. 

"d^" 

(^ ; Mode combined g^ (10'^) 90% UL low g^ (10~") 90% UL high g^ (10~^) 90% UL 

(N 

^ '. K°^j.+fj,- OAQtoil ±0.03 (2.2cr) 0.99 O.Oll^ 13 ±0.01 (O.Ocr) 0.31 0.251° il ±0.02 (1.9cr) 0.53 

■ K+^i+fj,~ 0-^'i-to,\t ±0.02 (S.Ocr) 0.62 O.l23to°on ±0.007 (l.Scr) 0.25 O.lietlfn ±0.010 (1.7cr) 0.24 

K°e+e- O.OSt^i^ ±0.01 (0.6cr) 0.36 -OM9to°o% ±0.004 (O.Ocr) 0.10 0.177l|5:Jgf ±0.011 (2.3cr) 0.38 

K+e+e- 0.5llo:n ±0.02 (5.8cr) - O.SOStHea ±0.013 (6.7cr) - 0.125to:o55 ±0.006 (2.5cr) 0.21 

><' 

(U ; K'+n+fi- 1.46to:75 ±0.12 (2.0a) 2.55 0.75i;;:4^ ±0.07 (1.70-) 1.48 0.78lJ!;4g ±0.06 {2.1a) 1.36 

P i' K*°jj.+fi'' 1.35to:3? ±0.10 (4.1cr) - 0.4lt°Jt ±0.04 (2.2cr) 0.83 0.62to:?9 ±0.05 (3.8o-) 0.97 

^ '. K'+e+e' 1.38+0-42 ±0.08 (3.7cr) 1.97 l-06+°;28 ±0-07 (4.7cr) - O.ig+oj? ±0.01 (0.9<t) 0.53 

K'^e+e- 0.86to:24 ±0-05 (4.2cr) - 0.20to:n ±0.01 (2.1cr) 0.42 0.3515513 ±0-02 (3.0o-) 0.60 



^ 



(N 

>, 

0^ , In addition to the above results, wc have also performed fits in extended q^ regions for K* di-clcctron final states, 

where a significant rate enhancement compared to that above q^ > 0.1 GeV^/c^ is expected. Table HIl shows the total 
and partial branching fractions for these modes for extended q^ regions. 

^^ ' TABLE IL Branching fractions for the extended q^ regions with statistical and systematic errors, respectively. 

o 



X 



Mode combined g^ (10"®) low g^ (10"'') 



C^ ■ B+-*K*+e+e- l.gotlj;^^ ± 0.11 1.32115;^^ ± 0.09 

B°^K*°e+e- 1.02tl52° ± 0.06 0.731 



We also perform fits combining various hadronic and di-lepton final states which are then averaged assuming isospin 
and/or lepton-flavor symmetry. Table Hill shows branching fractions for the combined modes including statistical and 
systematic errors, statistical significance and upper limits for measurements with < Aa statistical significance. 



TABLE III: Combined mode branching fractions by q^ region with statistical and systematic errors, respectively. The statistical 
significance of each measurement is given parenthetically. 90% confidence level upper limits are provided for measurements 
with < 4(T statistical significance. 



Mode 


combined q^ (10"^) 


90% UL 


low g^ (10-'^) 


90% UL 


high q^ (10-«) 


90% UL 


KH+r 


0.2lt°:\l ±0.02 (1.7(7) 


0.45 


-OMltoM7 ±0.004 (0.0(7) 


0.09 


0.202l°j;^? ±0.013 (3.0(7) 


0.36 


K+i+r 


0.476l[5:^^6 ±0.022 (6.5(t) 


- 


0.250to:o4? ±0.010 (6.7(7) 


- 


0.122to:o44 ±0.006 (3.0(7) 


0.19 


Kix+ti- 


OAltoll ±0.02 (4.70-) 


- 


0.096+0059 ±0.006 (2.1(7) 


0.21 


0.139tg;o6l ±0.008 (2.9(7) 


0.25 


Ke+e- 


0.3881^83 ±0.020 (5.6a) 


- 


0.217toMl ±0.010 (5.9(7) 


- 


0.132+°g|o ±0.007 (3.5(7) 


0.21 


Ki+r 


0.394155:°^^ ±0.020 (7.3cr) 


— 


0-18lto:o36 ±0.008 (6.1(7) 


— 


0.135t^;|;^? ±0.007 (4.5(7) 


— 


K*H+r 


imtoll ±0.07 (5.8(7) 




0-257J;;;:J^^ ±0.02 (2.9(7) 


0.44 


0.46l'^;l? ±0.03 (4.7(7) 




K*+£+r 


lAOtoi° ±0-09 (4.20-) 


- 


0.98l|^:2| ±0.06 (4.9(7) 


- 


0.34t'^;?J ±0.02 (1.8(7) 


0.63 


K*^^+fl- 


1.35l°:^^ ±0-10 (5.30-) 


- 


0.47l!^:ig ±0.04 (3.2(7) 


0.84 


0.65t|5;i? ±0.05 (4.7(7) 


- 


K*e+e~ 


0-99J:;5i? ±0.06 (5.6ct) 


- 


0A2to]_l ±0.03 (4.3(7) 


- 


0.30t°;l? ±0.02 (3.1(7) 


0.49 


K*e+r 


l.llt;;:^^ ±0.07 (7.7(7) 


— 


0.43lE;:lJ ±0.03 (5.3(7) 


— 


0.42l°/i° ±0.03 (5.3(7) 


— 



